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1. Historical Overview 

One could fill volumes (and many people have done so) writing on the history of super- 
fluidity and superconductivity. Our goals are much humbler: we only wish to connect the 
pairing found in neutron matter with recent cold fermionic atom experiments. Thus, our 
ever-shifting kaleidoscope (which relies mainly on Refs. [l]|2l[3l|4l) by construction passes 
over in silence discoveries that are major, but do not bear a direct relation to the physics 
of neutron-star crusts. The natural way to start such a brief historical overview is by men- 
tioning Heike Kamerlingh Onnes, who in 1908 liquefied ^He, at 4.2 K. Onnes mentions in 
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his 1913 Nobel lecture [5| that the behavior of ^He at an even colder temperature, 2.2 K, 
presents an interesting problem, which "could possibly be connected with the quantum the- 
ory" (note that this was said in 1913). Onnes had already realized that there was something 
peculiar about that temperature. Instead of staying on the subject of bosonic systems, he 
used "^He to cool down mercury, discovering superconductivity (a drop of the resistivity to 
non-measurable values) at 4.1 K in that system. 

Then, in 1925, when the statistics different particles follow was still unresolved, Ein- 
stein predicted the phenomenon now known as Bose-Einstein condensation (BEC), accord- 
ing to which below a certain temperature a finite fraction of the total number of particles 
occupies the lowest-energy single-particle state. In 1927, Keesom and Wolfke made the 
distinction between the two phases of ^He above and below the A transition point around 
2.2 K. In 1938 Kapitsa and, separately, Allen & Misener discovered superfluidity (a drop 
of the viscosity to essentially zero) using ^He. A year later, Kapitsa managed to get Lev 
Davidovich Landau released from prison, with the hope that he would be able to provide a 
theoretical explanation of this new effect. 

In an apparently unrelated context, Baade and Zwicky (in 1934, only two years after 
the discovery of the neutron) proposed the idea of high-density and small-radius neutron 
stars which would be much more gravitationally bound than ordinary stars and would be 
formed in supernova explosions. On the terrestrial front, Fritz London was the first person 
to suggest that superconductivity and superfluid flow in liquid helium were macroscopic 
quantum phenomena. Then in 1950 Ginzburg and Landau proposed what seemed to be a 
phenomenological theory by introducing a complex pseudo wave function ip as an order 
parameter in Landau's general theory of second-order phase transitions. In 1957 the epoch- 
making microscopic Bardeen-Cooper-Schrieffer (BCS) theory of superconductors was put 
forward, allowing one to calculate the pairing gap A. In it, the electron-phonon interaction 
causes an instability in the Fermi-sea with respect to the formation of what are now called 
Cooper pairs. Just two years later, Gorkov showed that the Ginzburg-Landau theory was 
a limiting form of the BCS theory around the critical temperature, in which ijj is directly 
proportional to the gap A. 

In the meantime, Bohr, Mottelson, and Pines had already (in 1958) made an analogy 
between the low-energy spectra of nuclei and of the electrons in a superconducting metal. 
In 1959 Migdal, ||6l in what was essentially a throw-away comment, noted that the "su- 
perfluidity of nuclear matter may lead to some interesting cosmological phenomena if stars 
exist which have neutron cores. A star of this type would be in a superfluid state with a tran- 
sition temperature corresponding to 1 MeV". Thus, the study of superfluidity in infinitely 
extended nuclear systems (such as neutron-star matter or nuclear matter) predates the 1967 
discovery of pulsars, which were identified as rapidly rotating magnetic neutron stars the 
next year. A few years after its formation, the relevant part of a neutron star is expected to 
have a temperature from 10^ K to 10^ K, i.e. a small fraction of the pairing gap. 

Enter unconventional superfluidity: in 1972 Lee, Osheroff, and Richardson observed 
superfluid phases of ^He (a fermionic atom) at a transition temperature of around 2.7 mK, 
three orders of magnitude lower than for ^He. Just weeks after the discovery of the super- 
fluid, Anthony Leggett extended previous work of his and identified the states seen in the 
new data. In a similarly untypical development, in 1986 Bednorz and Miiller discovered 
superconductivity at 35 K in a barium lanthanum copper oxide, initiating the era of what 
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are now known as high- Tc superconductors. Currently, the highest-Tc superconductor is 
mercury thallium barium calcium copper oxide at 138 K. In contradistinction to the case of 
■^He, a comprehensive theory of high-Tc superconductors is still being sought. 

In parallel developments, in 1980 Leggett showed that for fermions the limits of tightly 
bound molecules and long-range Cooper pairs are connected in a smooth crossover (BCS- 
BEC crossover). On the bosonic side, 1995 saw the first realization of Bose-Einstein con- 
densation in dilute vapors of alkali atoms. Most of the physics involved was governed by 
mean-field interactions, and could be well described by the Gross-Pitaevskii equation. As 
an experimental discovery it was very important in and for itself, providing a clear realiza- 
tion of the phenomenon that Einstein had theoretically predicted 70 years earlier But the 
methods involved in this research program on alkali atoms (magnetic trapping, evaporative 
cooling, sympathetic cooling, optical trapping, and Feshbach resonances) were then used to 
examine fermionic gases (in some cases, by the same research groups), the first example of 
quantum degeneracy in trapped Fermi gases being "^"K in 1999 at JILA. In November 2003, 
three groups reported the realization of Bose-Einstein condensates of molecules composed 
of fermionic atoms, and soon after that the observations were extended throughout the en- 
tire BEC-BCS crossover One of the high points (possibly the high point) in the ensuing 
explosion of research was the demonstration of fermionic superfluidity and phase coher- 
ence at MIT in April 2005. These phenomena, just like superfuidity in neutron matter, are 
characterized by strong interactions and cannot be dealt with satisfactorily in a mean-field 
theory framework. However, in contradistinction with neutron matter, cold atoms can be 
directly probed in the laboratory and can therefore help constrain nuclear theory. 

2. Mise-en-scene 

The materials already referred to as high-Tc superconductors are so called because they 
have the highest Tc's (in an absolute scale) currently known on earth. However, if we take 
a different view and speak instead in terms of relative temperatures, comparing the critical 
temperature with a temperature that sets the scale in each system (the Fermi temperature 
Tp), things look different. In Tabled] we have summarized some results for fermionic sys- 
tems. Neutron matter and ultracold atomic gases are separated by 17 orders of magnitude: 
however, in terms of their corresponding Fermi temperatures, their critical temperatures are 

Table 1. Typical values for the transition temperature Tc, and the ratio of the transition 
temperature to the Fermi temperature Tc/Tp, for various fermionic superfluids or super- 
conductors. 



Conventional superconductors 
^He 

High-Tc superconductors 
Neutron matter 
Atomic Fermi gases 



2.7 mK 
100 K 
IO^OK 
200 nK 



5 K 



5 -IQ-s 
5 -10-4 
10-2 



0.1 
0.2 
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very large, essentially of comparable magnitude. This is the reason both these systems can 
be viewed as "high-temperature superfluids", even though one occurs at very small and the 
other at very large temperature. Since the pairing gap A is roughly proportional to the crit- 
ical temperature Tc, both these systems also exhibit very strong pairing, the strongest ever 
observed before in nature or experimented with in the laboratory. As we shall see below, in 
cold fermionic atoms the particle-particle interactions can be tuned experimentally, allow- 
ing one to examine a whole range of new phenomena, thus also mimicking the setting of 
low-density neutron matter, which is beyond direct experimental reach. In order to make 
this more tangible, let us begin with Schwinger and Bethe's "shape-indendent formula", 
used to describe s-wave scattering: 



In this expression, k is the momentum and 5 is the s-wave phase shift. Knowing the latter 
gives us the scattering cross section, which is experimentally observable. This formula 
implies that low-energy scattering provides information about only these two parameters: 
the scattering length a and the effective range r^. Thus, the cross section is independent of 
the details of the potential. Let us now discuss some specifics of our two systems and their 
interparticle interactions. 

2.1. Ultracold Atomic Gases 

For ultracold fermionic atoms at low temperature, |[3]|4l|7ll3| the superfluid phase arises only 
in the presence of interactions. The first batch of fermionic experiments was done with an 
equal mixture of the two lowest hyperfine states, conventionally referred to as hyperfine 
states |1) and |2), of ^Li or confined optically in a dipole trap. 

The atom-atom interaction is quite complicated, though sophisticated models such as 
the Aziz potential for helium are often replaced by the simplified Lennard- Jones form. For 
cold alkali gases we have a van der Waals length rc ~ 50 — 100 A, but in the quantum 
degenerate regime the interparticle spacing is at least a few thousand A, i.e. the interparticle 
spacing is significantly larger than the characteristic length of the interaction. In a broad 
Feshbach resonance we have kprf. <C 1, where kp is the Fermi wave vector (connected 
to the number density via p = A;|./(37r^)), so at fixed density the effective-range can be 
taken to be very small, essentially zero (though it may be possible to use narrow and wide 
resonances in cold atoms to study the case of varying r^. experimentally, [9] and thus directly 
simulate neutron matter). 

What was critically new in the fermionic cold atom experiments of the last decade 
was the ability to vary the interaction strength (the scattering length a) across a resonance, 
through a regime known as "unitarity". This system is very interesting, since it is both 
dilute (i.e. the range of the interatomic potential is much smaller than the interparticle 
distance) and strongly interacting (i.e. the scattering length is much larger than the inter- 
particle distance). In this regime, all the length scales associated with interactions drop out 
of the problem and the system is expected to exhibit universal behavior, [10] i.e. a non- 
dependence on the specifics of the interatomic potential, which is why, in principle, any 
fermionic atom can be used to perform such experiments. The unitary Fermi gas has been 
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the subject of intensive theoretical study, both in the homogeneous case, and for the case of 
a trapped gas. 

Let us now move on to observable many-body properties that can be and have been 
measured in such experiments. Experiments using ^Li at Duke University 1 8 ] and at ENS 
lITl I have measured the ground-state energy of the system, essentially finding it to be in 
agreement with Quantum Monte Carlo predictions |[T2l [T3l [T4l [T5l [T6ll (very recent exper- 
imental results from MIT ifTTl and fermion-sign-free theoretical values ifTSl change this 
picture only quantitatively). The ground-state energy per particle is conventionally given in 
units of the energy of a free Fermi gas at the same density Epc = = kj, / (lOm) 

as E = ^EpG- Other experiments at MIT and Rice probed lithium gases with population 
imbalance (also called "polarized" gases). An MIT experiment fT9l established the phase 
diagram of a polarized gas, revealing spatial discontinuities in the spin polarization. This 
experiment was then used |20| to extract the pairing gap, which was found to be approxi- 
mately half of the Fermi energy Ep, in good agreement with QMC calculations |[T4|. (Sim- 
ilarly to the case of the ground-state energy, the gap is conventionally given in units of Ep 
as A = rjEp). The MIT group went on to use RF spectroscopy to independently determine 
the gap, finding it to be in agreement with the afore-mentioned calculation and extraction. 
ED 

2.2. Neutron-Star Matter 

The neutrons in the inner crust are expected to pair in the ^Sq channel, |[22ll23l but super- 
fluidity in this channel will probably exhibit gap closure before reaching the core. Let us 
note that even though, as we already pointed out, the temperature of the crust is from 10^ K 
to 10^ K the neutron gas is still superfluid because the critical temperature is expected to be 
larger, approximately 10^*^ K (see Table[Tl). The first-order approximation to this superfluid 
gas embedded in a lattice of nuclei is to consider infinite pure neutron matter and examine 
the effect of the ion lattice only at a later stage. For many decades, there has been a large 
spread of predictions of the ^Sq pairing gap even for this idealized pure system, which we 
discuss in a later section. This spread is mainly due to the differences in the many-body 
schemes employed: dilute neutron matter is a system that is very strongly paired, meaning 
that all perturbative approaches fail, and most diagrammatic techniques in general have no 
way of specifying why any specific class of diagrams should be dominant. 

Even so, dilute neutron matter is simple when compared to the general nuclear many- 
body problem. To see this, we turn turn our attention to the neutron-neutron interaction. 
In general it is complicated, necessitating sophisticated models (such as the Argonne vl8 
potential) to describe it. It has a strong dependence on quantum numbers such as total spin 
S and isospin T. In local potentials it has a strongly repulsive core at distances of 0.5 fm, 
as well as an electromagnetic (in the case of neutrons, magnetic-moment) term, and one- 
pion exchange at large distances and intermediate-range spin-dependent attraction which is 
dominated by two-pion exchange. However, as in the case of cold fermionic atoms, there are 
some simplifications: in very dilute neutron matter the density (or kp) goes to zero, making 
all the terms in the interaction apart from the s-wave channel drop off. In contradistinction 
to the case of cold atoms, the neutron-neutron effective range is not necessarily small in 
comparison to the intemeutron spacing: = 2.7 fm (with kpVe going from to 1.5). 
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Another critical distinction is that for neutrons in a neutron-star crust the scattering length 
is not at our disposal, being fixed, instead, at the value of a = —18.5 fm, larger than the 
interparticle spacing {kpa goes from to -10). 

Let us now move on to possible observational evidence of neutron-star matter pairing. 
Superfluidity in a neutron star is often used to explain its dynamical and thermal evolution. 
Starting with the former: the main example of this is given by pulsar post-glitch timing 
observations. For some pulsars, even though their angular frequency decreases slowly but 
steadily, there also occur sudden and rapid speed-ups of the frequency (called glitches), 
which are in their turn followed by a gradual decay of the increased frequency. The dom- 
inant phenomenological model holds that the speed-up is triggered by a "starquake" in the 
crust. The crust speed-up is rapidly communicated to the charged particles by the strong 
magnetic field. However, the response of the neutron superfluid to the speed-up is consid- 
erably slower and proceeds via electron scattering by the normal fluid cores of the vortex 
lines in the rotating superfluid. (This model was, astonishingly, proposed immediately after 
the first Vela pulsar ghtch.) lEU 

The other commonly referred to aspect in which pairing is important has to do with 
neutron-star cooling: observations of cooling quiescent neutron stars provide evidence for 
the existence of a neutron superfluid throughout the inner crust. ll25l This is so because the 
specific heat of a superfluid is exponentially suppresed. Furthermore, in the presence of a 
neutron ^Sq gap, the neutron-neutron bremsstrahlung reaction rate is also suppressed. Fi- 
nally, Cooper-pair breaking/formation neutrino emission processes that occur near the tran- 
sition temperature are also relevant to the cooling of neutron stars, though due to the small 
size of the crust this happens only during the first few hundred years of evolution, i.e. during 
the crust's thermal relaxation. Cooper-pair breaking/formation calculations have recently 
introduced vector current conservation, leading to the conclusion that the vector part of the 
contribution to this process is strongly suppressed, though the axial part is only slightly 
modified. These three factors (specific heat, bremsstrahlung, pair breaking/formation) have 
been incorporated in a recent minimal-cooling paradigm calculation, [26| which showed 
that different ^Sq pairing gap calculations lead to only slight quantitative differences in the 
cooling curves of young neutron stars. All these factors show the importance of superflu- 
idity in neutron matter, but do not probe the magnitude of the pairing gap. The situation 
in that regard may be different as regards a new mechanism that makes use of superfluid 
phonons. [27] Whether this mechanism is competitive to the heat conduction by electrons 
in magnetized neutron stars or not is a question that is directly correlated to the size of the 
gap- 
Finally, let us point out that, outside the observational realm, neutron matter computa- 
tions also hold significance in the context of traditional nuclear physics: equation of state 
results at densities close to the nuclear saturation density have been used for some time 
to constrain Skyrme and other density functional approaches to heavy nuclei, [28] while 
the density-dependence of the ^Sq gap in low-density neutron matter has recently also 
been used to constrain Skyrme-Hartree-Fock-Bogoliubov treatments in their description 
of neutron-rich nuclei ll29l . 
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3. Quantum Many-Body Methods 

In the series of works reviewed in the present text, both cold atoms and neutron matter 
were attacked using the same methodology: Quantum Monte Carlo many-body simulations. 
Such methods have been used for decades to calculate the energies of many systems, e.g. 
liquid helium. In such systems, however, QMC methods were unable to provide the pairing 
gaps, because of the vast difference in scale between the energy of the entire system (of the 
order of eV's per particle) and the pairing gap (of the order of meV). Thus, in a fortuitous 
turn of events, the same feature of strongly paired fermionic systems (the large pairing gap) 
that precludes the application of mean-field theories, is precisely the reason that allows 
many-body simulation techniques to be used. Before giving a few technical details on how 
Green's Function Monte Carlo works, we first go over the analytically exact expectations 
for quantum many-body theory, that we know any correct new approach should recover 
as a limiting case. After that we discuss BCS mean-field theory, which is not expected to 
be exact in the strongly paired regime, but even so is a useful benchmark with which to 
compare other theories. 



3.1. Weak coupling 

At extremely low densities i\kpa\ << 1) the effective coupling between two fermions 
is weak and matter properties can be calculated analytically. The ground-state energy of 
normal (i.e. non-superfluid) matter in this regime was calculated by Lee and Yang in 1957: 
AMI 

E . 10 , 4 .2 



_ 1 + — A:^a + —^ (11 -2 In 2) (A:i.a)^ (2) 

where Epc is the energy of a free Fermi gas at the same density as the interacting gas. 
While this expression ignores the contributions of superfluidity, these are exponentially 
small in (l/kpa). In the next section we compare these results to QMC calculations for 

Ikpa] > 1. 

The pairing gap at weak coupling is also known analytically. The mean-field BCS 
approach described below [Eq. |9l reduces in this limit to: 

However, as was shown in 1961 by Gorkov and Melik-Barkhudarov, ||3T1 the BCS result 
acquires a finite polarization correction even at weak coupling, yielding a reduced pairing 
gap: 

AO/, N 1 8 h^k^p / \ 

^^^^) = (4^?^^^H^J ■ 

Thus, the polarization corrections reduce the mean-field BCS result by a factor of 
l/(4e)^/^ ^ 0.45. Interestingly, if one treats the polarization effects at the level of sophis- 
tication used in the work of Gorkov and Melik-Barkhudarov, this factor changes with kpa 
1321, though there is no a priori reason to expect such an approach to be valid at stronger 
coupling (kpa of order 1 or more). Calculating the pairing gap in this region has been an 
onerous task, as can be seen from the multitude of publications devoted to this subject in 
the past few decades. |[22l[23|[33|[34l[Ml [36] [T7l[M|[t9l[40ll^ 
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3.2. BCS in the continuum and in a box 

As the coupling strength increases, we expect the BCS mean-field theory to become more 
accurate. In the BCS-BEC transition studied in cold atoms, the BCS theory goes correctly to 
the two-body bound state equation in the deep BEC regime. Though we do not expect BCS 
results to be quantitatively reliable in the regime of interest for neutron-star crusts, BCS 
theory provides a standard basis of comparison for our ab initio results and also allows us 
to analyze finite-size effects in the QMC simulations in a simple way. Within the BCS 
formalism the wave function is: 

iV')=nK+^^k4/_k4)io), (5) 

k 

where + 1;^ = 1. A variational minimization of the expectation value of the Hamiltonian 
for an average partice number (or density) leads to the gap equation: 



where the elementary quasi-particle excitations of the system have energy: 



(6) 



S(k) = ^e(k)2 + A(k)2 (7) 

and ^(k) = e(k) — ^, where the chemical potential is /i and e(k) = is the single- 
particle energy of a particle with momentum k. The chemical potential is found by solving 
the gap equation together with the equation that provides the average particle number: 



(8) 



When interested in the ^Sq gap for neutron matter, it is customary to perform partial-wave 
expansions of the potential and the gap functions, as well as an angle-average approxima- 
tion. Thus, Eq. ^ takes the form: 



oo 



A(A.) = -i/d^'fc-^A(^'), (9) 







where the potential matrix element is: 

oo 

v{k,k') = Jdrr'^jo{k'r)V{r)jo{kr). (10) 



Similarly, Eq. ^ becomes: 

oo 
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ni , 1 , 1 , 1 , 1 , 1 , 1 1 

2 4 6 8 10 12 

-kpa 

Figure 1. BCS pairing gap A divided with the Fermi energy Ep, versus the Fermi mo- 
mentum kp for cold atoms (solid line), neutrons with AVI 8 (dashed line), neutrons with a 
modified Poschl-Teller potential (dotted line) tuned to have the same scattering length and 
effective range as AVI 8. At low density all curves are identical for practical purposes. Also 
shown is the solution of the BCS problem in a periodic box using the modified Poschl-Teller 
potential for 66 neutrons (dash-dotted line). 



These equations are one dimensional, and thus simple to treat numerically. The density 
equation can be decoupled from the gap equation only when A/fi << 1; this is not the 
case for the density regime we are considering. 

In Ref. Us I Eq. Q was solved in tandem with Eq. (fTTI) for a modified Poschl-Teller 
(MPT) potential: 

v{r) = -vo -Y- — r , (12) 

m cosh (i/r) 

where and u are parameters which can be tuned so that this potential reproduces any 
desirable scattering length and effective range. The potential in Eq. ([T2l ) clearly has no 
repulsive core, making it amenable to a stiaightforward iterative solution. In this reference 
the potential was tuned so that the effective range is considerably shorter than the average 
interparticle spacing and at the same time the product kpa was varied from -1 to -10. The 
result is the solid line marked as "atoms" in Fig. [T] If kpa is taken to be infinite (the 
unitarity regime, unattainable in neutron star crusts) the BCS limiting value turns out to be 
0.69. 

In Ref. Il44l Eq. Q was then solved together with Eq. (ITTI ) using the same MPT 
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potential but the potential parameters were tuned so as to reproduce the neutron-neutron 
scattering length a « —18.5 fm and effective range rg ^ 2.7 fm, the result also being 
shown in Fig. [T] (dotted line). The same exercise was repeated for a more realistic micro- 
scopic potential, namely the ^Sq channel of the Argonne vl8 |45| potential that contains a 
strong short-range repulsion. This calculation is greatly simplified if one uses the method 
described in Ref. H6l . thereby transforming the problem into a quasilinear one; the result 
is the dashed line. For all the densities considered, the results of solving Eqs. ^ and ([TT]) 
with the appropriately tuned MPT potential and with the Argonne vl8 potential are virtually 
indistinguishable. 

The modified Poschl-Teller potential can also be used in a calculation for finite average 
particle number. This is interesting because the Quantum Monte Carlo simulations dis- 
cussed below are performed in a periodic simulation volume. This means that, in principle, 
there may exist considerable finite-size effects: since the simulation cannot be performed 
in the thermodynamic limit, the question arises of how close to it one is. Refs. |[T5l l44l |47l 
solved Eqs. Q and ^ for (N) from 20 to 200, in periodic boundary conditions in a cubic 
box of volume L^: 



There it was found that (A'') = 66 is very close to the thermodynamic limit. In Fig. [T] we 
show the results of solving the BCS gap equation Eq. Q in a periodic box along with the 
particle-number conserving Eq. ^ for {N) = 66 (dash-dotted line); as already noted, the 
discrepancy between that curve and the continuum result is small. Given this finite-size 
analysis, the Quantum Monte Carlo results shown in the rest of this review were performed 
using approximately 66 particles in periodic boundary conditions. 

3.3. Quantum Monte Carlo 

The Hamiltonian for cold atoms and neutron matter at low densities is: 



where N is the total number of particles. In the case of cold atoms the interaction only acts 
between opposite spin pairs and is of the form shown in Eq. (fT2l) . The neutron-neutron 
interaction is more complicated but still simple if one considers the AV4 formulation: HSl 



— inx,ny,nz) ■ 



(13) 




(14) 



V4{r) = Vc{r) + v^{r)(Ti ■ <T2 



(15) 



In the case of 5=0 (singlet) pairs this gives: 



vs{r) = Vc{r) - ?>v„{r) . 



(16) 



However, it also implies an interaction for S=\ (triplet) pairs: 



vp{r) = Vc{r) + v„{r) . 



(17) 



Ref. 1441 explicitly included such p-wave interactions in the same-spin pairs (the contribu- 
tion of which was small even at the highest density considered), and perturbatively corrected 
the 5 = 1, Ms = pairs to the correct p-wave interaction. 
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In these calculations it is customary to first employ a standard Variational Monte Carlo 
simulation. This approach uses Monte Carlo integration to minimize the expectation value 
of the Hamiltonian: 

^^^^^^^ - JdIl\^y{K)\^ - • ^^^^ 

thereby optimizing a variational wave function "^y 

At a second stage, the output of the Variational Monte Carlo calculation is used as input 
in a fixed-node Green's Function Monte Carlo simulation, which projects out the lowest- 
energy eigenstate "^q from a trial (variational) wave function "i>v- This it does by treating 
the Schrodinger equation as a diffusion equation in imaginary time r and evolving the 
variational wave function up to large r. The ground state is evaluated from: 

^-0 = exp[-{H - Et)t]'^v (19) 

evaluated as a branching random walk. The short-time propagator is usually taken as 

exp[-HAT] = exp[-yAr/2] exp[-rAr] exp[-V6T/2], (20) 

which is accurate to order (Ar)^. 

The fixed-node calculation gives a wave function "^q that is the lowest-energy state with 
the sames nodes (surface where ^' = 0) as the trial state ^y. The resulting energy Eq is 
an upper bound to the true ground-state energy. The variational wave function "^y has a 
Jastrow-BCS form (see below), and contains a variety of parameters, many of which affect 
the nodal surfaces. Since the fixed-node energy is an upper bound to the true ground state, 
these parameters can be optimized to give the best approximation to the ground-state wave 
function. In order to optimize these variational parameters, they were included as slowly 
diffusing coordinates in a preliminary GFMC calculation. The parameters evolve slowly 
in imaginary time, equilibrating around the lowest-energy state consistent with the chosen 
form of the trial wave function. [12 J 

The ground-state energy Eq can be obtained from: 

^ {^v\H\^o) ^ {^o\H\^o) 

In all these Quantum Monte Carlo superfluid simulations the trial wave function was 
taken to be of the Jastrow-BCS form with fixed particle number: 

= n fp(^^j) n fpi^i'f) n /(^uo^in ^(^if)] (22) 

and periodic boundary conditions. The primed (unprimed) indices correspond to spin-up 
(spin-down) neutrons. The pairing function (/)(r) is a sum over the momenta compatible 
with the periodic boundary conditions. In the BCS theory the pairing function is: 

Hr) = E ^e^^"- = E "ne*''"" > (23) 
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and can be parametrized with a short- and long-range part as in Ref. [12]: 

</.(r) = /3(r) + (24) 

n, I<Ic 

where / = n^+n^+nl using the parameters defined in Eq. ( [T3l ). The Jastrow part is usually 
taken from a lowest-order-constrained-variational method [49] calculation described by a 
Schrodinger-like equation: 

- —^^f{r) + v{r)f{r) = Xf{r) (25) 
m 

for the opposite-spin /(r) and by the corresponding equation for the same-spin fp{r). 
Since the /(r) and /p(r) are nodeless, they do not affect the final result apart from reducing 
the statistical error. The fixed-node approximation guarantees that the result for one set of 
pairing function parameters in Eq. (l24l) will be an upper bound to the true ground-state 
energy of the system. As mentioned above, the parameters are optimized in the full QMC 
calculation, providing the best possible nodal surface, in the sense of lowest fixed-node 
energy, with the given form of trial function. This upper-bound property is made use of in 
an attempt to get as close to the true ground-state energy as possible. 



4. Equation of state 

The T = equations of state for cold atoms and neutron matter of Ref. 1 15 1 are compared 
in Fig. 121 The horizontal axis is kpa, with the equivalent Fermi momentum kp for neutron 
matter shown along the top. The vertical axis is the ratio of the ground-state energy to the 
free Fermi gas energy {Epc) at the same density; it must go to one at very low densities 
and decrease as the density increases and the interactions become important. The curve 
at lower densities shows the analytical result by Lee and Yang we mentioned in an earlier 
section. The QMC results shown in this figure seem to agree with the trend implied by the 
Lee- Yang result (which has not been extended to even higher density because of its weak- 
coupling nature). The neutron matter and cold atom equations of state are very similar even 
for densities where the effective range is comparable to the interparticle spacing. Hence, 
our earlier hope that cold-atom experiments can tell us something rather directly about 
the neutron-matter equation of state is now more than just a hope. Near kpa = —10 
the energy per particle is not too far from QMC calculations and measurements of the 
ratio ^ of the unitary gas energy to Epc shown as an arrow on the right (it corresponds 
to kpa = oo); previous calculations give ^ k, 0.4 (however, as mentioned above, see also 
Refs. iHTl.fTSll). At larger densities the cold-atom and neutron matter results start to diverge 
from each other. This is due to two separate factors: i) the neutron finite effective range, and 
ii) the fact that the neutron results also incorporate an early attempt to correct the S = 1, 
Ms = pairs. When the density is very low, the s-wave contribution is dominant so cold 
atoms and neutron matter agree very well. At higher densities the energy is higher with the 
contribution of the p-wave interaction. 

We also compare our neutron matter QMC AV4 results from Ref. P4l for the ground- 
state energy to other calculations extending to larger Fermi momenta. Let us note that 
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Figure 2. Quantum Monte Carlo equation of state for cold atoms (squares) and neutron 
matter (circles). Also shown are the analytic expansion of the ground-state energy of a 
normal fluid (line) and the Quantum Monte Carlo result at unitarity (arrow). 



these are slightly different (larger) than the values shown in Fig. |2] because of i) includ- 
ing same-spin p-wave contributions, and ii) a more refined approach to dealing with the 
5 = 1, Ms = pairs. For the highest density examined, kpa = —10, this change is 
approximately 7%, while for kpa = —5 it is only 1%. Nonperturbative corrections at the 
highest density considered could reduce the difference between the s-wave interaction and 
AV4 results slightly. In Fig. [3]the results of Ref. [44] are compared to (approximate) vari- 
ational hypemetted-chain calculations by Friedman and Pandharipande ll50l . and another 
calculation by Akmal, Pandharipande, and Ravenhall (APR) [1511 . Also included are two 
Green's Function Monte Carlo results for 14 neutrons with more complete Hamiltonians 
||52| . a result following from a Brueckner-Bethe-Goldstone expansion [55], a difermion 
EFT result (shown are the error bands) [|53l . the latest Auxiliary-Field Diffusion Monte 
Carlo (AFDMC) result (discussed below) [HH . a Dirac-Brueckner-Hartree-Fock calcula- 
tion [54], a lattice chiral EFT method at next to leading order [56j| (see also Ref. [i571 ). and 
an approach that makes use of chiral N^LO three-nucleon forces. ll59l Of these, Refs. [ [5T1 . 
||4TI . and |[59l include a three-nucleon interaction, though at the densities considered, these 
are not expected to be significant. Qualitatively all of these results agree within 20%. 



14 



Alexandres Gezerlis 



kp [fm"'] 



0.2 0.4 0.6 0.8 1 



0.8 
0.6 
0.4 



' 1 ' 1 


1 


X 


Friedman-Pandharipande 






► 


APR 






▲ 


GFMC v6, v8' 


- 




• 


Baldo-Maieron 


^ Lee-Yang 






Schwenk-Pethick 




* 


XEFTNLO 






+ 


DBHF 






▼ 


Hebeler-Schwenk 


- ■ 




♦ 


AFDMC 






■ 


QMC AV4 


♦ 








■ ♦ 


X 








- 1 


X 


• 




* * ^ + ' f ^ 


1,1,1 


1 




1,1,1,1 



2.5 5 7.5 10 12.5 15 17.5 20 

- k„ a 

F 



Figure 3. Equation of state for neutron matter compared to various previous results. De- 
spite quantitative discrepancies, all calculations give essentially similar results. Our lowest 
density corresponds to kpa = —1. 



5. Pairing gap 

In Fig. |4]we plot the pairing gap as a function of kpa for both cold atoms and neutron 
matter, taken from Ref. 1 15 1. The Quantum Monte Carlo pairing gaps are arrived at from 
the ground-state energy, through the use of the odd-even staggering formula: 

A = E{N + 1) - ^ [EiN) + E{N + 2)] , (26) 

where N is an even number of particles. BCS calculations are shown as solid lines, and 
QMC results are shown as points with eiTor bars. QMC pairing gaps are shown from cal- 
culations of = 66 — 68 particles. For very weak coupling, —kpa « 1, the pairing gap 
is expected to be reduced from the BCS value by the polarization corrections calculated by 
Gorkov and Melik-Barkhudarov, A/A^cs = (l/4e)^/'^, as mentioned in an earlier section 
(Eq. Because of finite-size effects, it is difficult to calculate pairing gaps using QMC 
in the weak coupling regime. The QMC calculations at the lowest density, kpa = —1, are 
roughly consistent with this reduction from the BCS value. At slightly larger yet still small 
densities, where —kpa = 0{1) but kpr^ << 1 for neutron matter, one would expect the 
pairing gap to be similar for cold atoms and neutron matter. The results at kpa = —2.5, 
where kpr^, « 0.35, support this expectation. Beyond that density the effective range be- 
comes important and the QMC results are significantly reduced in relation to the cold atoms 
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Figure 4. Superfluid pairing gaps versus kpa for cold atoms {r^ ~ 0) and neutron matter 
(|re/a| f» 0.15). BCS (solid lines) and QMC results (points) are shown. 



where ~ 0. 

In cold atoms, the suppression from BCS is reduced as the density increases, with a 
smoothly growing fraction of the BCS results as we move from the BCS to the BEC regime. 
At unitarity the measured pairing gaps fT9l l20ll2r i are 0.45(0.05) of the Fermi energy, for 
a ratio A/Abcs ~ 0.65, in agreement with predictions by QMC methods. |[T2l [T4l [T5]l 
In neutron matter, though, the finite range of the potential reduces A/Ep as the density 
increases. We find a ratio A/Abcs that increases slightly from |fci?a| = 1 to 2.5, but 
thereafter remains roughly constant. 

In Fig. [5] we compaie the neutron matter results of Ref. ll44l to selected previous 
results: a CoiTclated-Basis Function calculation by Chen et al. |[33l . an extension of the 
polarization-potential model by Wambach et al. 041 . a medium-polarization calculation by 
Schulze et al. 1135 1 , a renormalization group calculation by Schwenk et al. 1361 . a Brueckner 
calculation by Cao et al. |39|, a determinantal lattice QMC approach Ii42l . and finally the 
newer Correlated-Basis Function calculation by Fabrocini et al. |[38i that was used as an 
input wave function in the two AFDMC calculations of 2005 and 2008. |[38ll4ll 

The results of our QMC calculations are much larger than most diagrammatic ll33l [34l 
[35l and renormalization group fW\ approaches. As these approaches assume a well-defined 
Fermi surface or calculate polarization corrections based on single-particle excitations it is 
not clear how well they can describe neutron matter in the strongly paired regime, or the 
similar pairing found in cold atoms. Ref. 1391 . which incorporates self-energy corrections 
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Figure 5. Superfluid pairing gap versus kpa for neutron matter compared to previous re- 
sults. 



and screening at the RPA level within Brueckner theory, appears to give results similar 
to ours. However, these values disagree with our lower-density results and, perhaps more 
importantly, at the lowest density reported the gap is larger than the mean-field BCS value 
(which is disallowed, as discussed in section 13.1.1 ). On a similar note, Refs. (34"] and |[36l 
make use of a weak-coupling formula to calculate the pairing gap, similarly to Eqs. (O and 
dUl. The prefactor they use is justified based on predictions in the theory of ^He. However, 
the concept itself of a Fermi surface is not well-defined in these strongly paired systems: in 
^He, in contrast to the present case, the gap is considerably smaller than the Fermi energy. 

Having discussed the relationship of our QMC results to analytical approaches, a few 
comments on the relationship between the different extant simulation results are in order. 
The first AFDMC results for pairing gaps in neutron matter were published in Ref. |38|. 
These initial calculations were for very small system sizes of N= 14-20 and yielded quite 
large pairing gaps. The AFDMC group then repeated their calculations for larger systems, 
141 1 reducing the gap and eliminating part of the difference with GFMC results, but not 
entirely. Detailed explanations are given in Ref. [44 J. The AFDMC approach includes a 
realistic two plus three nucleon interaction, and is therefore valuable for moderate to large 
densities. At very small densities we expect the s- and p-wave interactions to be sufficient. 
The trial wave function (and the constrained path approximation) for the AFDMC method 
are chosen from another approach. In the calculations of Refs. ll38ll4T1l the wave function 
was taken from a Correlated-Basis Function (CBF) approach. This wave function is unre- 
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alistic at very low density as it yields a pairing gap larger than the Fermi energy. The QMC 
AV4 results use a wave function that has been variationally optimized, but a much simpler 
Hamiltonian that should nevertheless be valid at the lowest densities. Ref. B4l demon- 
strates that QMC calculations of the gap can be quite sensitive to the trial function, for 
example using the CBF input wave function the QMC calculation reproduces the AFDMC 
result for the pairing gap. We expect the QMC results are more accurate at low densities, 
but at higher densities the full interaction treated with AFDMC is required, along with a 
careful study of the impact of the trial wave function. 

Finally, the QMC results seem to qualitatively agree (at least for the lowest densities 
considered) with a determinantal Quantum Monte Carlo lattice calculation [42J which, how- 
ever, only includes the s-wave component in the interaction (which is why they agree much 
better with the results of Ref. fl5 1 than those of Ref. [44]). This implies that a consensus is 
emerging, in that both the GFMC and the dQMC approaches find a gap that is suppressed 
with respect to the mean-field BCS result but is still a substantial fraction of the Fermi 
energy. 

6. Conclusions and Future Directions 

To conclude, we have calculated the equation of state and pairing gap of cold atoms and low- 
density neutron matter with the AV4 interaction from \kpa\ from 1 to 10. The calculated 
equation of state and pairing gap match smoothly with the known analytic results at low 
densities, and provide important constraints in the strong-coupling regime at large kpa. 
The low-density equation of state can help constrain Skyrme mean-field models of finite 
nuclei. The pairing gap for low-density neutron matter is relevant to Skyrme-Hartree-Fock- 
Bogoliubov calculations (29] of neutron-rich nuclei and to neutron-star physics, since it 
is expected to influence the behavior of the crust. [26] Moreover, a magnetic field in the 
neutron star crust would have to be approximately 10^^ G to overcome this gap and thus 
polarize neutron matter; such a value of the magnetic field is not implausible within the 
context of magnetars. The question of polarized neutron matter also has a long history ll60l 
and has recently been attacked using Quantum Monte Carlo. lIHTl l62l Similarly, the fact 
that the magnitude of the gap is not as small as previously expected implies that a new 
mechanism that makes use of superfluid phonons is competitive to the heat conduction by 
electrons in magnetized neutron stars. ll27l 

This line of Quantum Monte Carlo calculations, having first been applied to and ver- 
ified in cold atomic experiments, can also provide directions for future work in the field 
of nucleonic infinite matter. The simplest case is that of a two-component gas where the 
two populations are equal. ifTSl l44l The next step is to examine the ramifications of taking 
different populations for the two components: this is the case of spin-polarized low-density 
neutrons studied in Refs. lIHTl l62l . Cold-atom experiments have by now also addressed 
Efimov physics, in wliich three components are involved. In the nuclear context, adding 
a third species could provide further insight into the physics of neutron stars. If the third 
component particles were taken to be protons and, as in this paper, only a few of them 
were added, then it would be possible to study highly asymmetric nuclear matter. Another 
possible avenue of future research is related to optical lattice experiments with cold atoms: 
to first approximation these are equivalent to periodic external potentials. In the nuclear 
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case, an external potential would allow us to study the static response of neutron matter and 
would also facilitate the understanding of the impact on neutron pairing of the ion lattice 
that exists in a neutron star crust. Such microscopic results for the static response could pro- 
vide further constraints on energy-density functionals used to describe the crust of neutron 
stars. 
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